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SUMMARY 

Numerical calculations of the 2-D steady incompressible driven cavity flow are presented. The Navier- 
Stokes equations in streamfunction and vorticity formulation are solved numerically using a fine 
uniform grid mesh of 601 x 601. The steady driven cavity solutions are computed for Re < 21,000 
with a maximum absolute residuals of the governing equations that were less than 10 -10 . A new 
quaternary vortex at the bottom left corner and a new tertiary vortex at the top left corner of the 
cavity are observed in the flow field as the Reynolds number increases. Detailed results are presented 
and comparisons are made with benchmark solutions found in the literature. 
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1. INTRODUCTION 

Numerical methods for 2-D steady incompressible Navier-Stokes (N-S) equations are often 
tested for code validation, on a very well known benchmark problem; the lid-driven cavity flow. 
Due to the simplicity of the cavity geometry, applying a numerical method on this flow problem 
in terms of coding is quite easy and straight forward. Despite its simple geometry, the driven 
cavity flow retains a rich fluid flow physics manifested by multiple counter rotating recirculating 
regions on the corners of the cavity depending on the Reynolds number. In the literature, it is 
possible to find different numerical approaches which have been applied to the driven cavity 
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flow problem [1-34]. Though this flow problem has been numerically studied extensively, still 
there are some points which are not agreed upon. For example; 1) an interesting point among 
many studies is that different numerical solutions of cavity flow yield about the same results for 
Re < 1,000 however, start to deviate from each other for larger Re, 2) also another interesting 
point is that while some studies predict a periodic flow at a high Reynolds number, some 
others present steady solutions for even a higher Reynolds number. The objective of this work 
is then to investigate these two points and present accurate very fine grid numerical solutions 
of steady 2-D driven cavity flow. 

Among the numerous studies found in the literature, we will give a brief survey about some 
of the significant studies that uses different types of numerical methods on the driven cavity 
flow. In doing so the emphasis will be given on three points; on the numerical method used, 
on the spatial order of the numerical solution and on the largest Reynolds number achieved. 

Recently, Barragy & Carey [2] have used a p-type finite element scheme on a 257 x 257 
strongly graded and refined element mesh. They have obtained a highly accurate (Ah 8 order) 
solutions for steady cavity flow solutions up to Reynolds numbers of Re=12,500. Although 
Barragy & Carey [2] have presented qualitative solution for Re=16,000 they concluded that 
their solution for Re=16,000 was under-resolved and needed a greater mesh size. 

Botella & Peyret [6] have used a Chebyshev collocation method for the solution of the lid- 
driven cavity flow. They have used a subtraction method of the leading terms of the asymptotic 
expansion of the solution of the N-S equations in the vicinity of the corners, where velocity 
is discontinuous, and obtained a highly accurate spectral solutions for the cavity flow with a 
maximum of grid mesh of N = 160 (polynomial degree) for Reynolds numbers Re < 9,000. 
They stated that their numerical solutions exhibit a periodic behavior beyond this Re. 

Schreiber & Keller [28] have introduced an efficient numerical technique for steady viscous 
incompressible flows. The non-linear differential equations are solved by a sequence of Newton 
and chord iterations. The linear systems associated with the Newton iteration is solved by 
L [/-factorization with partial pivoting. Applying repeated Richardson extrapolation using the 
solutions obtained on different grid mesh sizes (maximum being 180 x 180), they have presented 
high-order accurate (Ah 8 order in theory) solutions for Reynolds numbers Re < 10,000. 

Benjamin & Denny [5] have used a method of relaxing the algebraic equations by means of 
the ADI method, with a non-uniform iteration parameter. They have solved the cavity flow 
for Re < 10,000 with three different grid mesh sizes (maximum being 101 x 101) and used a 
Ah n extrapolation to obtain the values when Ah — > 0. 

Wright & Gaskell [36] have applied the Block Implicit Multigrid Method (BIMM) to the 
SMART and QUICK discretizations. They have presented cavity flow results obtained on a 
1024 x 1024 grid mesh for Re < 1,000. 

Nishida & Satofuka [24] have presented a new higher order method for simulation of the 
driven cavity flow. They have discretized the spatial derivatives of the N-S equations using a 
modified differential quadrature (MDQ) method. They have integrated the resulting system 
of ODEs in time with fourth order Runge-Kutta-Gill (RKG) scheme. With this they have 
presented spatially Ah 10 order accurate solutions with grid size of 129 x 129 for Re < 3,200. 

Liao [20], and Liao & Zhu [21], have used a higher order streamfunction-vorticity boundary 
element method (BEM) formulation for the solution of N-S equations. With this they have 
presented solutions up to Re=10,000 with a grid mesh of 257 x 257. 

Hou et. al. [17] have used Lattice Boltzmann Method for simulation of the cavity flow. They 
have used 256 x 256 grid points and presented solutions up to Reynolds numbers of Re=7,500. 
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Goyon [13] have solved the streamfunction and vorticity equations using Incremental 
Unknowns. They have presented steady solutions for Re < 7,500 on a maximum grid size 
of 256 x 256. 

Rubin & Khosla [27] have used the strongly implicit numerical method with 2x2 coupled 
streamfunction-vorticity form of the N-S equations. They have obtained solutions for Reynolds 
numbers in the range of Re < 3,000 with a grid mesh of 17 x 17. 

Ghia et. al. [12] later have applied a multi-grid strategy to the coupled strongly implicit 
method developed by Rubin & Khosla [27]. They have presented solutions for Reynolds 
numbers as high as Re=10,000 with meshes consisting of as many as 257 x 257 grid points. 

Gupta [15] has used a fourth (Ah 4 ) order compact scheme for the numerical solution of 
the driven cavity flow. He has used a 9 point (3x3) stencil in which the streamfunction and 
vorticity equations are approximated to fourth order accuracy. He has used point-SOR type of 
iteration and presented steady cavity flow solutions for Re < 2,000 with a maximum of 41 x 41 
grid mesh. 

Li et. al. [19] have used a fourth (Ah 4 ) order compact scheme which had a faster convergence 
than that of Gupta [15]. They have solved the cavity flow with a grid size of 129 x 129 for Re 
< 7,500. 

Bruneau & Jouron, [7] have solved the N-S equations in primitive variables using a full 
multigrid-full approximation storage (FMG-FAS) method. With a grid size 256 x 256, they 
have obtained steady solutions for Re < 5,000. 

Grigoriev & Dargush [14] have presented a BEM solution with improved penalty function 
technique using hexagonal subregions and they have discretized the integral equation for each 
subregion as in FEM. They have used a non-uniform mesh of 5040 quadrilateral cells. With 
this they were able to solve driven cavity flow up to Re=5,000. 

Aydin & Fenner [1] have used BEM formulation with using central and upwind finite 
difference scheme for the convective terms. They have stated that their formulation lost its 
reliability for Reynolds numbers greater than 1,000. 

To the authors best knowledge, in all these studies among with other numerous papers found 
in the literature, the maximum Reynolds number achieved for the 2-D steady incompressible 
flow in a lid-driven cavity is Re=12,500 and is reported by Barragy & Carey [2]. Although 
Barragy & Carey [2] have presented solutions for Re=16,000, their solution for this Reynolds 
number display oscillations related to boundary layer resolution issues due to a coarse 
mesh. Nallasamy & Prasad [22] have presented steady solutions for Reynolds numbers up 
to Re < 50,000, however, their solutions are believed to be inaccurate as a result of excessive 
numerical dissipation caused by their first order upwind difference scheme. 

A very brief discussion on computational as well as experimental studies on the lid-driven 
cavity flow can be found in Shankar & Deshpande [29]. 

Many factors affect the accuracy of a numerical solution, such as, the number of grids in the 
computational mesh (Ah), and the spatial discretization order of the finite difference equations, 
and also the boundary conditions used in the solution. 

It is an obvious statement that as the number of grids in a mesh is increased (smaller Ah) 
a numerical solution gets more accurate. In this study the effect of number of grid points 
in a mesh on the accuracy of the numerical solution of driven cavity flow is investigated, 
especially as the Reynolds number increases. For this, the governing N-S equations are solved 
on progressively increasing number of grid points (from 128 x 128 to 600 x 600) and the solutions 
are compared with the highly accurate benchmark solutions found in the literature. 
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Order of spatial discretization is another factor on the accuracy of a numerical solution. 
Finite element and spectral methods provide high spatial accuracy. This study also investigates 
the effect of spacial accuracy on the numerical solution of driven cavity flow. For this, following 
[28] , repeated Richardson extrapolation is used and then the extrapolated results are compared 
with the benchmark solutions. 

Another factor on the accuracy of a numerical solution is the type and the order of the 
numerical boundary conditions used in computation. For driven cavity flow, this subject was 
briefly discussed by Weinan & Jian-Guo [34], Spotz [30], Napolitano et. al. [23] and also by 
Gupta & Manohar [16]. 

Until Barragy & Carey [2], many studies have presented steady solutions of driven cavity 
flow for Re=10,000 ([28], [5], [21] and [12]). They ([2]) have presented solutions for Re=12,500. 
In all these studies that have presented high Reynolds number numerical solutions, the largest 
grid size used was 256 x 256. In this study, the grid size will be increased up to 600 x 600 and 
the effect of the grid size on the largest computable Reynolds number solution of the driven 
cavity flow will be investigated. 

This paper presents accurate, very fine grid (600 x 600) numerical solutions of 2-D steady 
incompressible flow in a lid-driven cavity for Reynolds numbers up to Re — 21, 000. A detailed 
comparisons of our results with mainly the studies mentioned above as well as with other 
studies not mentioned here, will be done. 



2. NUMERICAL METHOD 

For two-dimensional and axi-symmetric flows it is convenient to use the streamfunction (ip) 
and vorticity (w) formulation of the Navier-Stokes equations. In non-dimensional form, they 
are given as 

d 2 ip d 2 tjj 

lh 2 + dy 2= ~ UJ U 

1 / d 2 uj d 2 cu \ dip duj dip dm 
Re \ dx 2 By 2 J dy dx dx dy 

where, Re is the Reynolds number, and x and y are the Cartesian coordinates. The N-S 
equations are non-linear, and therefore need to be solved in an iterative manner. In order to 
have an iterative numerical algorithm, pseudo time derivatives are assigned to these equations 
(1 and 2) such as 

dtp d 2 ip d 2 ip 

duj 1 d 2 uj 1 d 2 Lo dip duj dip dco 
dt Re dx 2 Re dy 2 dx dy dy dx 

Since we are seeking a steady solution, the order of the pseudo time derivatives will not 
affect the final solution. Therefore, an implicit Euler time step is used for these pseudo time 
derivatives, which is first order (At) accurate. Also for the non-linear terms in the vorticity 
equation, a first order (At) accurate approximation is used and the governing equations (3 
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and 4) are written as the following 

dx 2 " dy 2 



l/} n+i_ At ^ -At^ =tp n + Atuj n (5) 



u n+1 - At — — - At— — ^ -At 9 ^-"— ^ At 9 ^"— n+1 -co n (6) 
Re dx 2 Re dy 2 dx dy dy dx 

In operator notation the above equations look like the following. 

1 - A, S - A 'w) 5 '" + ' - + A '"" (7) 

*-"is£-"iS&-M5)"£ + "(2)"s) »" + ' ^ < 8 > 

The above equations are in fully implicit form where each equation requires the solution of 
a large banded matrix which is not computationally efficient. Instead of solving the equations 
(7 and 8) in fully implicit form, more efficiently, we spatially factorize the equations, such as 

' 1 - A, £)( 1 - i '^)*" +, = *" + A "'" (9) 
+ Af (I)" I) (' - "k $ - « it) I) - ™> 

Here, we split the single LHS operators into two operators, where each only contains derivatives 
in one direction. The only difference between equations (7 and 8) and (9 and 10) are the At 2 
terms produced by the factorization. The advantage of this process is that each equation now 
require the solution of a tridiagonal system, which is numerically more efficient. 

As the solution converges to the steady state solution, the At 2 terms due to the factorization 
will not become zero and will remain in the solution. To illustrate this, equation (9) can be 
written in explicit form as 

a2,/.n+l p*2 in+1 ai.i.n+1 

r +1 - - ^-t^- + At2 V^ = r + Atu> n (ii) 

dx z dy A dx z dy z 

As the solution converges to a steady state, \|/ n + 1 becomes equal to \P" and they cancel each 
other from the above equation and at convergence, equation (11) therefore appears as the 
following 

d 2 ij} d 2 ip ^ d 4 ip 

dx 2 dy 2 dx 2 dy 2 

In equation (12), there is no guarantee that the RHS, which is the consequence of the 
factorization, will be small. Only a very small time step (At) will ensure that the RHS of 
equation (12) will be close to zero. This, however, will slow down the convergence. In addition, 
the product may not be small since g^2g y 2 can be O(^) or higher in many flow situations. 

In order to overcome this problem, we add a similar term at the previous time level to the 
RHS of equation (11) to give the following finite difference form 

ffi.i.n+l pp. .I.n+l r) i ih n+1 fl A th n 



+ ld+"< = - a <I^ (12) 
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As a result, at steady state, ip n+1 converges to ip n and 9 a J^ dy2 converges to ^Jq^ , so that 
they cancel each out, and the final solution converges to the correct physical form 

d 2 ib d 2 ib , 

d^ + W + ^° (14) 

We would like to emphasize that the added At 2 term on the RHS is not an artificial numerical 
diffusion term. Such terms are often used to stabilize the numerical schemes so that convergence 
can be achieved. Here, the added At 2 term on the RHS is not meant to stabilize the numerical 
scheme. The only reason for this approach is to cancel out the terms introduced by factorization 
so that the equations are the correct physical representation at the steady state. 

One of the most popular spatially factorized schemes is undoubtedly the Beam & Warming 
[4] method. In that method, the equations are formulated in delta (A) form. The main 
advantage of a delta formulation in a steady problem is that there will be no second order 
(At 2 ) terms due to the factorization in the solution at the steady state. However in our 
formulation the second order (At 2 ) terms due to factorization will remain even at the steady 
state. Adding a second order (At 2 ) term to the RHS of the equations to cancel out the terms 
produced by factorization, easily overcomes this problem. One of the main differences between 
the presented formulation and a delta formulation is in the numerical treatment of the vorticity 
equation. In a delta formulation of the vorticity equation, there appears At order convection 
and dissipation terms on the RHS (explicit side) of the numerical equations. In the presented 
formulation, on the RHS of the equation there are no At order terms. Instead cross derivative 
terms ( dx 2 d 2 , gj.g , d ® d 2 , and d ® g ) appear on the RHS. Note that these terms are At 2 
order. Since At is small, the contribution of these At 2 terms are even smaller. Our extensive 
numerical tests showed that this approach has better numerical stability characteristics and 
is more effective especially at high Reynolds numbers compared to a delta formulation of the 
steady streamfunction and vorticity equations. 

Applying the approach of removing second order terms due to factorization on the RHS, to 
both equations (9 and 10), the final form of the numerical formulation we used for the solution 
of streamfunction and vorticity equation becomes 

1-A^) (l-A^) r+1 = r + A^ 

+ { At &){ At w) r (15) 
1 - + At (T\ I) d -Atl*- At W 9 



Re dx 2 \dy J dx ) \ Re dy 2 \ dx J dy 

+ { At Redx- 2 - At {i) dx){ At Redy- 2+At {i) dy) " (W) 

The solution methodology for the two equations involves a two-level updating. For the 
streamfunction equation, the variable / is introduced such that 

1 - Af w) r+1 = f (17) 
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and 

1 " ^ f = r+ At ^ + { At &) ( A V) < 18 > 

In equation (18) / is the only unknown. First, this equation is solved for / at each grid point. 
Following this, the streamfunction (ip) is advanced into the new time level using equation (17). 
In a similar fashion for the vorticity equation, the variable g is introduced such that 



and 



Re dx 2 \9y J d. 



As with /, g is determined at every grid point using equation (20), then vorticity (lj) is 
advanced into the next time level using equation (19). 

We note that, in this numerical approach the streamfunction and vorticity equations are 
solved separately. Each equation is advanced into a next time level by solving two tridiagonal 
systems, which allows the use of very large grid sizes easily. The method proved to be very 
effective on flow problems that require high accuracy on very fine grid meshes (Erturk & Corke 
[9] Erturk, Haddad & Corke [10]). 



3. RESULTS AND DISCUSSION 

The boundary conditions and a schematics of the vortices generated in a driven cavity flow 
are shown in Figure 1. In this figure, the abbreviations BR, BL and TL refer to bottom 
right, bottom left and top left corners of the cavity, respectively. The number following these 
abbreviations refer to the vortices that appear in the flow, which are numbered according to 
size. 

We have used the well-known Thorn's formula [32] for the wall boundary condition, such 
that 



^o = T^-VI (21) 



-2^i 2U 
AhF ~ Ah 

where subscript refers to points on the wall and 1 refers to points adjacent to the wall, Ah 
refers to grid spacing and U refers to the velocity of the wall with being equal to 1 on the 
moving wall and on the stationary walls. 

We note an important fact that, it is well understood ([34], [30], [18], [23]) that, even though 
Thorn's method is locally first order accurate, the global solution obtained using Thorn's 
method preserve second order accuracy. Therefore in this study, since three point second order 
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central difference is used inside the cavity and Thorn's method is used at the wall boundary 
conditions, the presented solutions are second order accurate 0(Ah 2 ). 

During our computations we monitored the residual of the steady streamfunction and 
vorticity equations (1 and 2) as a measure of the convergence to the steady state solution, 
where the residual of each equation is given as 

H 4> ~ + —9. + ( 22 ) 



Ax 2 Ay 2 



1 io n+ } ■ - 2lj u+1 + uj n i} ■ 1 uj n+1 , - 2uj n+1 + uj n+ h 



O _ <^±2± bl b JT±2l _|_ 

u Re Ax 2 Re Ay 



2 



1 „/,"+! , ,n+l , ,n+l „/,«+! „/,"+! , , 



2Ay 2Ax " " 2Aa; 2Ay ' 1 ' 

The magnitude of these residuals are an indication of the degree to which the solution has 
converged to steady state. In the limit these residuals would be zero. In our computations, for 
all Reynolds numbers, we considered that convergence was achieved when for each equation (22 
and 23) the maximum of the absolute residual in the computational domain (max(|_R^|) and 
max(|i? w |)) was less than 10 _10 . Such a low value was chosen to ensure the accuracy of the 
solution. At these residual levels, the maximum absolute difference in streamfunction value 
between two time steps, (max(|V>™ +1 - V> n |)), 

was in the order of 10 16 and for vorticity, 
(max(|o)" +1 — uj n \)), it was in the order of 10~ 14 . And also at these convergence levels, 
between two time steps the maximum absolute normalized difference in streamfunction, 
(maxd ^" 4 ^" I)), and in vorticity, (max(| |)), were in the order of 10" 13 , and 10" 12 
respectively. 

We have started our computations by solving the driven cavity flow from R= 1,000 to 
Re=5,000 on a 129 x 129 grid mesh. With using this many number of grids, we could not 
get a steady solution for Re=7,500. The first natural conclusion was that the cavity flow 
may not be steady at Re=7,500, and therefore steady solution may not be computable. Even 
though we have used a pseudo-time iteration, the time evolution of the flow parameters, either 
vorticity or streamfunction values at certain locations, were periodic suggesting that the flow 
indeed may be periodic, backing up this conclusion. Besides, in the literature many studies 
claim that the driven cavity flow becomes unstable between a Reynolds number of 7,000 and 
8,000, for example, Peng et. al. [25], Fortin et. al. [11] and Poliashenko & Aidun [26] predict 
unstable flow beyond Re of 7402, 7998,5 and 7763 respectively. However, on the other hand 
there were many studies ([2], [5], [12], [21], and [28]) that have presented solutions for a higher 
Reynolds number, Re=10,000. This was contradictory. We then have tried to solve the same 
case, Re=7,500, with a larger grid size of 257 x 257. This time we were able to obtain a 
steady solution. In fact, with 257 x 257 number of grids, we carry our steady computations 
up to Re=12,500. Barragy & Carey [2] have also presented steady solutions for Re=12,500. 
Interesting point was, when we tried to increase Re furthermore with the same grid mesh, 
we have again obtained a periodic behavior. Once again the natural conclusion was that the 
steady solutions beyond this Re may not be computable. Having concluded this, we could 
not help but ask the question what happens if we increase the number of grids furthermore. 
In their study Barragy & Carey [2] have presented a steady solution for Re=16,000. Their 
solution for this Reynolds number displayed oscillations. They concluded that this was due to 
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a coarse mesh. This conclusion encouraged us to use larger grid size than 257 x 257 in order to 
obtain solutions for Reynolds numbers greater than 12,500. At this point, since the focus was 
on the number of grids, we realized that among the studies that presented steady solutions at 
very high Reynolds numbers ([2], [5], [12], [21], and [28]), the maximum number of grids used 
was 257 x 257. Larger grid sizes haven't been used at these high Reynolds numbers. This was 
also encouraging us to use more grid points. We, then, decided to increase the number of grids 
and tried to solve for Re=15,000 with 401 x 401 number of grids. This time we were able to 
obtain a steady solution. This fact suggests that in order to obtain a numerical solution for 
2-D steady incompressible driven cavity flow for Re > 12,500, a grid mesh with more than 
257 x 257 grid points is needed. With using 401 x 401 grid points, we continued our steady 
computations up to Re=21,000. Beyond this Re our computations again displayed a periodic 
behavior. This time the first natural thing came to mind was to increase the number of grids. 
We have tried to use 513 x 513 grids, however again we could not get a steady solution beyond 
Re=21,000. Thinking that the increase in number of grids may not be enough, we then again 
increased the number of grids and have used 601 x 601 grids. The situation was the same, and 
the maximum Reynolds number that we can obtain a steady solution with using 601 x 601 
grids was Re=21,000. We did not try to increase the number of grids furthermore since the 
computations became time consuming. Whether or not steady computations of driven cavity 
flow are possible beyond Re=21,000 with grid numbers larger than 601 x 601, is still an open 
question. 

One of the reasons, why the steady solutions of the driven cavity flow at very high Reynolds 
numbers become computable when finer grids are used, may be the fact that as the number 
of grids used increases, Ah gets smaller, then the cell Reynolds number or so called Peclet 
number defined as Re c = decreases. This improves the numerical stability characteristics 
of the numerical scheme used ([31], [34]), and allows high cavity Reynolds numbered solutions 
computable. Another reason may be that fact that finer grids would resolve the corner vortices 
better. This would, then, help decrease any numerical oscillations that might occur at the 
corners of the cavity during iterations. 

We note that all the figures and tables present the solution of the finest grid size of 601 x 601 
unless otherwise stated. 

The accuracy of a finite difference solution is set by the mesh size, and by the spatial order 
of the finite difference equations and the boundary approximations. At low Reynolds numbers 
(Re=l,000), different numerical method solutions found in the literature agree with each other. 
However, the solutions at higher Reynolds numbers (Re> 1,000) have noticeable discrepancies. 
We believe it is mainly due to different spatial orders and different grid mesh sizes and different 
boundary conditions used in different studies. 

Table 1 tabulates the minimum streamfunction value, the vorticity value at the center of 
the primary vortex and also the center location of the primary vortex for Re=l,000 along with 
similar results found in the literature with the most significant ones are underlined. In Table 1 
among the most significant (underlined) results, Schreiber & Keller [28] have used Richardson 
extrapolation in order to achieve high spatial accuracy. For Re=l,000 their extrapolated 
solutions are Ah 6 order accurate and these solutions are obtained by using repeated Richardson 
extrapolation on three different mesh size solutions. 

Following the same approach, we have solved the cavity flow on three different grid mesh 
(401 x 401, 513 x 513, and 601 x 601) separately for every Reynolds number (Re=l,000 to 
Re=21,000). Since the numerical solution is second order accurate 0(Ah 2 ), the computed 
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values of streamfunction and vorticity have an asymptotic error expansion of the form 



We can use Richardson extrapolation in order to obtain high-order accurate approximations 
(Schreiber & Keller [28]). We have used repeated Richardson extrapolations using the three 
different mesh size solutions. The extrapolated values of the streamfunction and the vorticity at 
the center of the primary vortex are given in Table 2 and Table 3 respectively. The extrapolated 
results tabulated in Table 2 and Table 3 are, in theory, Ah 6 order accurate. 

Looking back to Table 1, for Re=l,000, our extrapolated results are in very good agrement 
with the extrapolated results of Schreiber & Keller [28] . In fact our extrapolated streamfunction 
value at the primary vortex differs only 0.013 % and the vorticity value differs only 0.095 % 
from the extrapolated results of Schreiber & Keller [28] . 

With spectral methods, very high spatial accuracy can be obtained with a relatively smaller 
number of grid points. For this reason, the tabulated results from Botella & Peyret [6] are 
believed to be very accurate. Our results are also in very good agreement such that our 
extrapolated results differs 0.016 % in the streamfunction value and 0.098 % in the vorticity 
value from the results of Botella & Peyret [6] . 

Barragy & Carey [2] have only tabulated the location of the vortex centers and the 
streamfunction values at these locations yet they have not presented any vorticity data. 
Nevertheless, their solutions based on p-type finite element scheme also have very high spatial 
accuracy (Ah s order accurate). They have presented the maximum Reynolds number solutions 
(Re=12,500) for the cavity flow, to the best of our knowledge, found in the literature. We will 
show that the agreement between our results and results of Barragy & Carey [2] is excellent 
through the whole Reynolds number range they considered (1,000 < Re < 12,500). Comparing 
solutions for Re=l,000 in Table 1, the difference is 0.022 % in the streamfunction value. 

Wright & Gaskell [36] have used a very fine grid mesh with 1024 x 1024 grid points in their 
solutions. With this many grids their solution should be quite accurate. Their streamfunction 
and vorticity value are very close to our solutions with a difference of 0.113 % and 0.114 % 
respectively. 

Nishida & Satofuka [24] have presented higher order solutions in which for Re= 1,000 their 
solutions are Ah 8 order accurate. Comparing with our results we find a very good agreement 
such that our extrapolated results differs 0.040 % in streamfunction value and 0.136 % in 
vorticity value. 

Grigoriev & Dargush [14] have solved the governing equations with a boundary element 
method (BEM). They have only tabulated the streamfunction values and the location of the 
centers of the primary and secondary vortices for driven cavity flow for Reynolds numbers up 
to 5,000. Their solution for Re=l,000 is quite accurate. 

Benjamin & Denny [5] have used a Ah n extrapolation to obtain the values when Ah — > 
using the three different grid mesh size solutions. Their extrapolated results are close to other 
underlined data in Table 1. It is quite remarkable that even a simple extrapolation produces 
such good results. 

Hou et. al. [17] have simulated the cavity flow by the Lattice Boltzmann Method. They have 
presented solutions for Re < 7,500 obtained by using 256 x 256 lattice. However their results 



ip(xi, y 3 ) = ipij + cxAh 2 + c 2 Ah i + ... 




(24) 
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have significant oscillations in the two upper corners of the cavity. Especially their vorticity 
contours have wiggles aligned with lattice directions, for Re > 1,000. 

Li et. al. [19] have used a compact scheme and presented fourth order accurate (Ah 4 ) results 
for a 129 x 129 grid mesh. Ghia et. al. [12] have used a second order method (Ah 2 ) with a mesh 
size of 129 x 129. These results, together with that of Goyon [13] and Liao & Zhu [21], might be 
considered as somewhat under- resolved. Our calculations showed that for a second order (Ah 2 ) 
spatial accuracy, even a 401 x 401 grid mesh solutions can be considered as under- resolved for 
Re=l,000. 

The results of Bruneau & Jouron, [7] and Vanka [33] appear over-diffusive due to upwind 
differencing they have used. 

From all these comparisons we can conclude that even for Re=l,000 higher order 
approximations together with the use of fine grids are necessary for accuracy. 

Figures 2 through 6 show streamfunction contours of the finest (601 x 601) grid solutions, 
for various Reynolds numbers. These figures exhibit the formation of the counter-rotating 
secondary vortices which appear as the Reynolds number increases. 

Table 4 compares our computed streamfunction and vorticity values at the center of the 
primary vortex with the results found in the literature for various Reynolds numbers. The 
agreement between our Richardson extrapolated results (Ah 6 order accurate) with results 
of Barragy & Carey [2] based on p-type finite element scheme (Ah s order accurate) is quite 
remarkable. One interesting result is that the streamfunction value at the center of the primary 
vortex decreases at first as the Reynolds number increases until Re=10,000, however then starts 
to increase as the Reynolds number increases further. This behavior was also documented by 
Barragy & Carey [2]. 

In terms of quantitative results, Gupta & Manohar [16], stated that for the numerical 
solutions of driven cavity flow, the value of the streamfunction (ip) at the primary vortex 
center appears to be a reliable indicator of accuracy, where as the location of the vortex center 
(a; and y) is limited by the mesh size, and the value of the vorticity (u>) at the vortex center 
is sensitive to the accuracy of the wall boundary conditions. Nevertheless, the locations of the 
primary and secondary vortices, as well as the values of the streamfunction (?/>) and vorticity 
(uj) at these locations are presented in Table 5. These are comparable with the solutions of 
Barragy & Carey [2] and Ghia et. al. [12]. It is evident that as the Reynolds number increases, 
the center of the primary vortex moves towards the geometric center of the cavity. Its location 
is almost invariant for Re > 17,500. Our computations indicate the appearance of a quaternary 
vortex at the bottom left corner (BL3) at Reynolds number of 10,000. We would like to note 
that, this quaternary vortex (BL3) did not appear in our solutions for Re=10,000 with grid 
sizes less than 513 x 513. This would suggest that, especially at high Reynolds numbers, in 
order to resolve the small vortices appear at the corners properly, very fine grids have to be 
used. Moreover, our solutions indicate another tertiary vortex at the top left side of the cavity 
(TL2) at Re=12,500. These vortices (BL3 and TL2) were also observed by Barragy & Carey 
[2]- 

Figures 7 and 8 present the u- velocity profiles along a vertical line and the v- velocity profiles 
along a horizontal line passing through the geometric center of the cavity respectively and also 
Tables 6 and 7 tabulates u- and v-velocity values at certain locations respectively, for future 
references. These profiles are in good agrement with that of Ghia et. al. [12] shown by symbols 
in Figures 7 and 8. 

We have solved the incompressible flow in a driven cavity numerically. A very good 
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mathematical check on the accuracy of the numerical solution would be to check the continuity 
of the fluid, as suggested by Aydin & Fenner [1], although this is done very rarely in driven 
cavity flow papers. We have integrated the u-velocity and v-velocity profiles, considered in 
Figures 7 and 8, along a vertical line and horizontal line passing through the geometric center 
of the cavity, in order to obtain the net volumetric flow rate, Q, through these sections. Since 
the flow is incompressible, the net volumetric flow rate passing through these sections should 
be equal to zero, Q = 0. The u-velocity and v-velocity profiles at every Reynolds number 
considered are integrated using Simpson's rule, then the obtained volumetric flow rate values 
are divided by a characteristic flow rate, Q c , which is the horizontal rate that would occur in 
the absence of the side walls (Plane Couette flow), to help quantify the errors, as also suggested 
by Aydin & Fenner [1] . In an integration process the numerical errors will add up. Nevertheless, 

f 1 udyl I f 1 vdx\ 

the volumetric flow rate values (Qi = -^^q — - and Q 2 = -^^q — ) tabulated in Table 8 are 
close to zero such that, even the largest values of Qi = 0.000000705 and Q 2 = 0.000002424 in 
Table 8 can be considered as Q\ « Q 2 ~ 0. This mathematical check on the conservation of 
the continuity shows that our numerical solution is indeed very accurate. 

Distinctively than any other paper, Botella & Peyret [6] have tabulated highly accurate 
Chebyshev collocation spectral vorticity data from inside the cavity, along a vertical line and 
along a horizontal line passing through the geometric center of the cavity, for Re=l,000. 
Comparing our solutions with theirs ([6]) in Figure 9 and 10, we find that the agreement with 
Botella & Peyret [6] is remarkable, with the maximum difference between two solutions being 
0.18%. 

Between the whole range of Re=l,000 and the maximum Re found in the literature, 
Re=12,500, our computed results agreed well with the published results. Since there is no 
previously published results (to the best of our knowledge) of steady driven cavity flow beyond 
Re > 12,500 in the literature, we decided to use Batchelor's [3] model in order to demonstrate 
the accuracy of our numerical solutions at higher Reynolds numbers. As seen in Figures 7 and 
8, the u- and v-velocity profiles change almost linearly in the core of the primary vortex as 
Reynolds number increases. This would indicate that in this region the vorticity is uniform. 
As Re increases thin boundary layers are developed along the solid walls and the core fluid 
moves as a solid body with a uniform vorticity, in the manner suggested by Batchelor [3]. At 
high Reynolds numbers the vorticity at the core of the eddy is approximately constant and 
the flow in the core is governed by 

dhjj_ + cfjj_ = c 

dx 2 dy 2 

where C is the constant vorticity value. With streamfunction value being ip = on the 
boundaries, inside the domain (x,y)=([-0.5,0.5],[-0.5,0.5]) the following expression is a solution 
to equation (25) 



( 2 1 °° \ 
h \^ A n cosh(2n — l)irx cos(2ra — l)iry ] 
2 8 / 



(26) 



where 



4f_l)n+l 1 

A n = J ±-L^( C os H n--)nr (27) 
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The mean square law proposed by Batchelor [3] relates the core vorticity value with the 
boundary velocities. From using the relation 

> [(vg 2 + (vg 2 ]ds = ^ |[/| 2 d S = i (28) 

where U is the boundary velocity, we can have 

11 16 ^ i 

& ~ 6 ~ ^ (2n-l) 4 (cosh(n- i)7r) 2 

+2 J ' ^ +7r^(-l)"(2n- l)^„cosh(2n- 1)ttx^ d.T (29) 

From the above expression, the numerical value of C (vorticity value at the core) can be 
found approximately as 1.8859, in accordance with the theoretical core vorticity value of 1.886 
at infinite Re, calculated analytically by Burggraf [8] in his application of Batchelor's model 
[3], which consists of an inviscid core with uniform vorticity, coupled to boundary layer flows 
at the solid surface. 

In Figure 11 we have plotted the variation of the vorticity value at the center of the 
primary vortex, tabulated in Table 2, as a function of Reynolds number. In the same figure, 
the theoretical value of 1.8859 is shown with the dotted line. This figure clearly shows that 
our computations, in fact, in strong agrement with the mathematical theory (Batchelor [3], 
Wood [35], and Burggraf [8]) such that computed values of the vorticity at the primary vortex 
asymptote to the theoretical infinite Re vorticity value as the Reynolds number increases. 



4. CONCLUSIONS 

Fine grid numerical solutions of the steady driven cavity flow for Reynolds numbers up 
to Re=21,000 have been presented. The steady 2-D incompressible Navier-Stokes equations 
in streamfunction and vorticity formulation are solved computationally using the numerical 
method described. The streamfunction and vorticity equations are solved separately. For each 
equation, the numerical formulation requires the solution of two tridiagonal systems, which 
allows the use of large grid meshes easily. With this we have used a fine grid mesh of 601 x 601. 
At all Reynolds numbers, the numerical solutions converged to maximum absolute residuals 
of the governing equations that were less than 10 -10 . 

Our computations indicate that fine grid mesh is necessary in order to obtain a steady 
solution and also resolve the vortices appear at the corners of the cavity, as the Reynolds 
number increases. Several unique features of the driven cavity flow have been presented such 
that the appearance of a quaternary vortex at the bottom left corner (BL3) at Reynolds number 
of 10,000 on fine grid mesh, and a tertiary vortex at the top left corner (TL2) at Reynolds 
number of 12,500. Detailed results were presented and they were compared with results found 
in the literature. The presented results were found to agree well with the published numerical 
solutions and with the theory as well. 
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Figure 1 



Schematic view of driven cavity flow 



Re=1,000 

Uniform Grid 601x601 




Figure 2. Streamline contours of primary and secondary vortices, Re=l,000 
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Figure 3. Streamline contours of primary and secondary vortices, Re=5,000 



Re=10,000 
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Figure 4. Streamline contours of primary and secondary vortices, Re=l 0,000 
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Figure 5. Streamline contours of primary and secondary vortices, Re=l 5,000 
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Figure 6. Streamline contours of primary and secondary vortices, Re=20,000 
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Figure 7. Computed u-velocity profiles profiles along a vertical line passing through 
the geometric center of the cavity at various Reynolds numbers 




Figure 8. Computed v- velocity profiles along a horizontal line passing through 
the geometric center of the cavity at various Reynolds numbers 




Figure 9. Computed vorticity values along a vertical line passing through 
the geometric center of the cavity, Re=l,000 




Figure 10. Computed vorticity values along a horizontal line passing through 
the geometric center of the cavity, Re=l,000 




Figure 11. Strength of vorticity at the primary vortex as a function of Reynolds number. 
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Table 1. Comparison of the properties of the primary vortex; the maximum streamfunction value, 
the vorticity value and the location of the center, for Re=1000 
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Table 2. Richardson extrapolation for streamfunction values 

at the primary vortex 



Re Grid u 0{Ah 2 ) u 0{Ah A ) u 0(Ah 6 ) 

1000 401x401 -2.062761 

513x513 -2.064765 -2.067904 

601x601 -2.065530 -2.067579 -2.067213 

2500 401x401 -1.961660 

513x513 -1.967278 -1.976078 

601x601 -1.969675 -1.976096 -1.976117 

5000 401x401 -1.909448 

513x513 -1.921431 -1.940201 

601x601 -1.926601 -1.940451 -1.940732 

7500 401x401 -1.878523 

513x513 -1.896895 -1.925673 

601x601 -1.904883 -1.926282 -1.926969 

10000 401x401 -1.853444 

513x513 -1.878187 -1.916945 

601x601 -1.888987 -1.917919 -1.919018 

12500 401x401 -1.830890 

513x513 -1.862010 -1.910757 

601x601 -1.875618 -1.912072 -1.913557 

15000 401x401 -1.809697 

513x513 -1.847199 -1.905943 

601x601 -1.863618 -1.907602 -1.909476 

17500 401x401 -1.789348 

513x513 -1.833212 -1.901921 

601x601 -1.852447 -1.903975 -1.906294 

20000 401x401 -1.769594 

513x513 -1.819761 -1.898343 

601x601 -1.841814 -1.900891 -1.903767 

21000 401x401 -1.761828 

513x513 -1.814492 -1.896986 

601x601 -1.837672 -1.899768 -1.902909 



Table 3. Richardson extrapolation for vorticity values 
at the primary vortex 





Re 
2500 


r AAA 

5000 


7500 


10000 


12500 


Present 






i> 

Ul 


-0.121470 
-1.976117 


-0.122233 
-1.940732 


-0.122386 
-1.926969 


-0.122390 
-1.919018 


-0.122326 
-1.913557 


Barragy & Carey 


Ref. 


[2] 


i> 

Ul 


-0.1214621 


-0.122219 


-0.1223803 


-0.1223930 


-0.1223584 


Schreiber & Keller 


Ref. 


[26] 


i> 

Ul 


- 


- 


- 


-0.12292 
-1.9263 


- 


Li et. al. 


Ref. 


[18] 


i> 

Ul 


- 


-0.120359 
-1.92430 


-0.1193791 
-1.91950 




- 


Benjamin & Denny 


Ref. 


[5] 

L J 


i> 

Ul 


- 


- 


- 


-0.1212 


- 


Ghia et. al. 


Ref. 


[12] 


1> 

Ul 


- 


-0.118966 
-1.86016 


-0.119976 
-1.87987 


-0.119731 
-1.88082 


- 


Liao & Zhu 


Ref. 


[20] 


1> 

Ul 


- 


-0.1186 
-1.9375 


-0.1201 
-1.9630 


-0.1201 
-1.9426 


- 


Hou et. al. 


Ref. 


[17] 


1> 

Ul 




-0.1214 


-0.1217 






Grigoriev & Dargush 


Ref. 


[14] 


1> 

Ul 




-0.12209 









Table 4. Comparison of the minimum streamfunction value, the vorticity value 
at the center of the primary vortex for various Reynolds numbers 



Re 

1000 2500 5000 7500 10000 12500 15000 17500 20000 21000 



Primary i/> 
Vortex uj 

(x , y) 


-0.118781 
-2.065530 
(0.5300 , 0.5650) 


-0.121035 
-1.969675 
(0.5200 , 0.5433) 


-0.121289 
-1.926601 
(0.5150 , 0.5350) 


-0.120924 
-1.904883 
(0.5133 , 0.5317) 


-0.120403 
-1.888987 
(0.5117 , 0.5300) 


-0.119831 
-1.875618 
(0.5117 , 0.5283) 


-0.119240 
-1.863618 
(0.5100 , 0.5283) 


-0.118641 
-1.852447 
(0.5100 , 0.5267) 


-0.118039 
-1.841814 
(0.5100 , 0.5267) 


-0.117797 
-1.837672 
(0.5100 , 0.5267) 


i> 

BRl u 

(x , y) 


0.17281B-02 
1.115505 
(0.8633 , 0.1117) 


0.26561E-02 
1.929029 
(0.8350 , 0.0917) 


0.30604E-02 
2.724481 
(0.8050 , 0.0733) 


0.32102E-02 
3.215956 
(0.7900 , 0.0650) 


0.31758E-02 
3.696991 
(0.7767 , 0.0600) 


0.30839E-02 
4.260493 
(0.7617 , 0.0550) 


0.29897E-02 
4.837680 
(0.7483 , 0.0500) 


0.28938E-02 
5.372747 
(0.7367 , 0.0467) 


0.27967E-02 
5.882504 
(0.7267 , 0.0450) 


0.27598E-02 
6.080586 
(0.7217 , 0.0433) 


t/> 

BLl lj 

(x , y) 


0.23261B-03 
0.353473 
(0.0833 , 0.0783) 


0.92541E-03 
0.966506 
(0.0850 , 0.1100) 


0.13639E-02 
1.502628 
(0.0733 , 0.1367) 


0.15171E-02 
1.838156 
(0.0650 , 0.1517) 


0.15930E-02 
2.150847 
(0.0583 , 0.1633) 


0.16355E-02 
2.318733 
(0.0550 , 0.1683) 


0.16484E-02 
2.441871 
(0.0533 , 0.1717) 


0.16331E-02 
2.594288 
(0.0517 , 0.1750) 


0.16032E-02 
2.823121 
(0.0483 , 0.1817) 


0.15897E-02 
2.877385 
(0.0483 , 0.1833) 


V> 

BR2 ui 

(x , y) 


-0.54962E-07 
-0.77076E-02 
(0.9917 , 0.0067) 


-0.12967E-06 
-0.15959E-01 
(0.9900 , 0.0100) 


-0.14010E-05 
-0.34108E-01 
(0.9783 , 0.0183) 


-0.30998E-04 
-0.158392 
(0.9517 , 0.0417) 


-0.13397E-03 
-0.302615 
(0.9350 , 0.0667) 


-0.24687E-03 
-0.397507 
(0.9283 , 0.0817) 


-0.33148E-03 
-0.453339 
(0.9267 , 0.0883) 


-0.39539E-03 
-0.512637 
(0.9283 , 0.0967) 


-0.44959E-03 
-0.544236 
(0.9300 , 0.1033) 


-0.47028E-03 
-0.574871 
(0.9317 , 0.1083) 


V> 

BL2 o> 

(x , y) 


-0.84221E-08 
-0.29802E-02 
(0.0050 , 0.0050) 


-0.32321E-07 
-0.93800E-02 
(0.0067 , 0.0067) 


-0.72080E-07 
-0.12252E-01 
(0.0083 , 0.0083) 


-0.20202E-06 
-0.16862E-01 
(0.0117 , 0.0117) 


-0.10151E-05 
-0.30981E-01 
(0.0167 , 0.0200) 


-0.57409E-05 
-0.68557E-01 
(0.0250 , 0.0317) 


-0.20172E-04 
-0.134120 
(0.0367 , 0.0417) 


-0.44887E-04 
-0.192907 
(0.0483 , 0.0483) 


-0. 74543 E-04 
-0.232067 
(0.0567 , 0.0533) 


-0.86415E-04 
-0.250751 
(0.0600 , 0.0550) 


V> 

TLl a; 

(x , y) 




0.34455E-03 
1.321269 
(0.0433 , 0.8900) 


0.14416E-02 
2.062775 
(0.0633 , 0.9100) 


0.21119E-02 
2.201598 
(0.0667 , 0.9133) 


0.25870E-02 
2.247599 
(0.0717 , 0.9117) 


0.29394E-02 
2.296811 
(0.0750 , 0.9117) 


0.32110E-02 
2.329198 
(0.0767 , 0.9117) 


0.34298E-02 
2.358627 
(0.0800 , 0.9133) 


0.36105E-02 
2.380582 
(0.0817 , 0.9133) 


0.36749E-02 
2.387858 
(0.0817 , 0.9133) 


1> 

BR3 ui 

(x , y) 


- 




0.21562E-09 
0.70828E-03 
(0.9983 , 0.0017) 


0.14409E-08 
0.22180E-02 
(0.9967 , 0.0033) 


0.51934E-08 
0.25744E-02 
(0.9967 , 0.0050) 


0.99032E-08 
0.38075E-02 
(0.9950 , 0.0050) 


0.13947E-07 
0.33588E-02 
(0.9950 , 0.0050) 


0.19492E-07 
0.43303E-02 
(0.9950 , 0.0067) 


0.29610E-07 
0.51574E-02 
(0.9933 , 0.0067) 


0.34929E-07 
0.63746E-02 
(0.9933 , 0.0083) 


BL3 ui 

(x , y) 


- 


- 






0.14382E-09 
0.78229E-03 
(0.0017 , 0.0017) 


0.50930E-09 
0.53920E-03 
(0.0017 , 0.0017) 


0.11280E-08 
0.10824E-02 
(0.0017 , 0.0033) 


0.22898E-08 
0.22499E-02 
(0.0033 , 0.0033) 


0.34584E-08 
0.18414E-02 
(0.0033 , 0.0033) 


0.37846E-08 
0.17128E-02 
(0.0033 , 0.0033) 


V> 

TL2 a; 

(x , y) 












-0.12752E-05 
-0.215534 
(0.0067 , 0.8317) 


-0.13394E-04 
-0.515658 
(0.0150 , 0.8267) 


-0.36218E-04 
-0.711567 
(0.0200 , 0.8233) 


-0.63084E-04 
-0.858460 
(0.0233 , 0.8200) 


-0.74273E-04 
-0.933479 
(0.0250 , 0.8200) 



Table 5. Properties of primary and secondary vortices; streamfunction and vorticity values, (x,y) locations 



y 


Re 
1000 


2500 


5000 


7500 


10000 


12500 


15000 


17500 


20000 


21000 


1.000 


1.0000 


1 


0000 


1 


0000 


1 


0000 


1 


0000 


1 


0000 


1.0000 


1 


0000 


1 


0000 


1 


0000 


0.990 


0.8486 





7704 





6866 





6300 





5891 





5587 


0.5358 





5183 





5048 





5003 


0.980 


0.7065 





5924 





5159 





4907 





4837 





4833 


0.4850 





4871 





4889 





4895 


0.970 


0.5917 





4971 





4749 





4817 





4891 





4941 


0.4969 





4982 





4985 





4983 


u.you 


n no 





4607 





4739 





4860 





4917 





4937 


U.4yd i 





4925 





4906 





4897 


0.950 


0.4582 





4506 





4738 





4824 





4843 





4833 


0.4811 





4784 





4754 





4742 


0.940 


0.4276 





4470 





4683 





4723 





4711 





4684 


0.4653 





4622 





4592 





4580 


0.930 


0.4101 





4424 





4582 





4585 





4556 





4523 


0.4492 





4463 





4436 





4425 


0.920 


0.3993 





4353 





4452 





4431 





4398 





4366 


0.4338 





4312 





4287 





4277 


0.910 


0.3913 





4256 





4307 





4275 





4243 





4216 


0.4190 





4166 





4142 





4132 


0.900 


0.3838 





4141 





4155 





4123 





4095 





4070 


0.4047 





4024 





4001 





3992 


0.500 


-0.0620 


-0 


0403 


-0 


0319 


-0 


0287 


-0 


0268 


-0 


0256 


-0.0247 


-0 


0240 


-0 


0234 


-0 


0232 


0.200 


-0.3756 


-0 


3228 


-0 


3100 


-0 


3038 


-0 


2998 


-0 


2967 


-0.2942 


-0 


2920 


-0 


2899 


-0 


2892 


0.180 


-0.3869 


-0 


3439 


-0 


3285 


-0 


3222 


-0 


3179 


-0 


3146 


-0.3119 


-0 


3096 


-0 


3074 


-0 


3066 


0.160 


-0.3854 


-0 


3688 


-0 


3467 


-0 


3406 


-0 


3361 


-0 


3326 


-0.3297 


-0 


3271 


-0 


3248 


-0 


3239 


0.140 


-0.3690 


-0 


3965 


-0 


3652 


-0 


3587 


-0 


3543 


-0 


3506 


-0.3474 


-0 


3446 


-0 


3422 


-0 


3412 


0.120 


-0.3381 


-0 


4200 


-0 


3876 


-0 


3766 


-0 


3721 


-0 


3685 


-0.3652 


-0 


3622 


-0 


3595 


-0 


3585 


0.100 


-0.2960 


-0 


4250 


-0 


4168 


-0 


3978 


-0 


3899 


-0 


3859 


-0.3827 


-0 


3797 


-0 


3769 


-0 


3758 


0.080 


-0.2472 


-0 


3979 


-0 


4419 


-0 


4284 


-0 


4142 


-0 


4054 


-0.4001 


-0 


3965 


-0 


3936 


-0 


3925 


0.060 


-0.1951 


-0 


3372 


-0 


4272 


-0 


4491 


-0 


4469 


-0 


4380 


-0.4286 


-0 


4206 


-0 


4143 


-0 


4121 


0.040 


-0.1392 


-0 


2547 


-0 


3480 


-0 


3980 


-0 


4259 


-0 


4407 


-0.4474 


-0 


4490 


-0 


4475 


-0 


4463 


0.020 


-0.0757 


-0 


1517 


-0 


2223 


-0 


2633 


-0 


2907 


-0 


3113 


-0.3278 


-0 


3412 


-0 


3523 


-0 


3562 


0.000 


0.0000 





0000 





0000 





0000 





0000 





0000 


0.0000 





0000 





0000 





0000 



Table 6. Tabulated u- velocity profiles along a vertical line passing through 
the geometric center of the cavity at various Reynolds numbers 



X 



Re 

1000 2500 5000 7500 10000 12500 15000 17500 20000 21000 



1.000 





0000 





0000 





0000 





0000 





0000 





0000 





0000 





0000 





0000 





0000 


0.985 


-0 


0973 


-0 


1675 


-0 


2441 


-0 


2991 


-0 


3419 


-0 


3762 


-0 


4041 


-0 


4269 


-0 


4457 


-0 


4522 


0.970 


-0 


2173 


-0 


3725 


-0 


5019 


-0 


5550 


-0 


5712 


-0 


5694 


-0 


5593 


-0 


5460 


-0 


5321 


-0 


5266 


0.955 


-0 


3400 


-0 


5192 


-0 


5700 


-0 


5434 


-0 


5124 


-0 


4899 


-0 


4754 


-0 


4664 


-0 


4605 


-0 


4588 


0.940 


-0 


4417 


-0 


5603 


-0 


5139 


-0 


4748 


-0 


4592 


-0 


4534 


-0 


4505 


-0 


4482 


-0 


4459 


-0 


4449 


0.925 


-0 


5052 


-0 


5268 


-0 


4595 


-0 


4443 


-0 


4411 


-0 


4388 


-0 


4361 


-0 


4331 


-0 


4300 


-0 


4287 


0.910 


-0 


5263 


-0 


4741 


-0 


4318 


-0 


4283 


-0 


4256 


-0 


4221 


-0 


4186 


-0 


4153 


-0 


4122 


-0 


4110 


0.895 


-0 


5132 


-0 


4321 


-0 


4147 


-0 


4118 


-0 


4078 


-0 


4040 


-0 


4005 


-0 


3975 


-0 


3946 


-0 


3936 


0.880 


-0 


4803 


-0 


4042 


-0 


3982 


-0 


3938 


-0 


3895 


-0 


3859 


-0 


3828 


-0 


3800 


-0 


3774 


-0 


3764 


0.865 


-0 


4407 


-0 


3843 


-0 


3806 


-0 


3755 


-0 


3715 


-0 


3682 


-0 


3654 


-0 


3627 


-0 


3603 


-0 


3593 


0.850 


-0 


4028 


-0 


3671 


-0 


3624 


-0 


3574 


-0 


3538 


-0 


3508 


-0 


3481 


-0 


3457 


-0 


3434 


-0 


3425 


0.500 





0258 





0160 





0117 





0099 





0088 





0080 





0074 





0069 





0065 





0063 


0.150 





3756 





3918 





3699 





3616 





3562 





3519 





3483 





3452 





3423 





3413 


0.135 





3705 





4078 





3878 





3779 





3722 





3678 





3641 





3608 





3579 





3567 


0.120 





3605 





4187 





4070 





3950 





3885 





3840 





3801 





3767 





3736 





3725 


0.105 





3460 





4217 





4260 





4137 





4056 





4004 





3964 





3929 





3897 





3885 


0.090 





3273 





4142 





4403 





4337 





4247 





4180 





4132 





4093 





4060 





4048 


0.075 





3041 





3950 





4426 





4495 





4449 





4383 





4323 





4273 





4232 





4218 


0.060 





2746 





3649 





4258 





4494 





4566 





4563 





4529 





4484 





4438 





4420 


0.045 





2349 





3238 





3868 





4210 





4409 





4522 





4580 





4602 





4601 





4596 


0.030 





1792 





2633 





3263 





3608 





3844 





4018 





4152 





4254 





4332 





4357 


0.015 





1019 





1607 





2160 





2509 





2756 





2940 





3083 





3197 





3290 





3323 


0.000 





0000 





0000 





0000 





0000 





0000 





0000 





0000 





0000 





0000 





0000 



Table 7. Tabulated v-velocity profiles along a horizontal line passing through 
the geometric center of the cavity at various Reynolds numbers 



Re 


n _ iJo 1 ^ 

Vl ~~ Qc 


_ \SaVdx\ 


1000 


0.000000045 


0.000000134 


2500 


0.000000046 


0.000000344 


5000 


0.000000067 


0.000000693 


7500 


0.000000090 


0.000001027 


10000 


0.000000114 


0.000001341 


12500 


0.000000143 


0.000001630 


15000 


0.000000163 


0.000001894 


17500 


0.000000259 


0.000002133 


20000 


0.000000557 


0.000002346 


21000 


0.000000705 


0.000002424 



Table 8. Volumetric flow rates through a vertical line, Qi, and a horizontal 
line, Q2, passing through the geometric center of the cavity 



